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ABSTRACT 

The widely existing shallow decay phase of the X-ray afterglows of gamma- 
ray bursts (GRBs) is generally accepted to be due to long-lasting energy injec- 
tion. The outflows carrying the injecting energy, based on the component that is 
dominative in energy, fall into two possible types: baryon-dominated and lepton- 
dominated ones. The former type of outflow could be ejecta that is ejected during 
the prompt phase of a GRB and consists of a series of baryonic shells with a distri- 
bution of Lorentz factors, and the latter type could be an electron-positron-pair 
wind that is driven by the post-burst central engine. We here provide a unified 
description for the dynamics of fireballs based on these two types of energy injec- 
tion, and calculate the corresponding high-energy photon emission by consider- 
ing synchrotron radiation and inverse Compton scattering (including synchrotron 
self-Compton and combined inverse-Compton) of electrons. We find that, in the 
two energy-injection models, there is a plateau (even a hump) in high-energy 
light curves during the X-ray shallow decay phase. In particular, a considerable 
fraction of the injecting energy in the lepton-dominated model can be shared 
by the long-lasting reverse shock since it is relativistic. Furthermore, almost 
all of the energy of the reverse shock is carried by leptons, and thus the inverse- 
Compton emission is enhanced dramatically. Therefore, this model predicts more 
significant high-energy afterglow emission than the baryon-dominated model. We 
argue that these observational signatures would be used to discriminate between 
different energy-injection models in the upcoming Gamma-ray Large Area Space 
Telescope (GLAST) era. 

Subject headings: gamma rays: bursts — radiation mechanism: nonthermal 
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1. Introduction 

As discovered by Swift, there is a shallow decay phase (temporal indices a ~ [0, —0.8]) 
from post-bm^st several tens of seconds to several hours (even days) occurring in the X- 
ray afterglow light curves of a significant fraction of gamma-ray bursts (GRBs) (Nousek et 
al. 2006; O'Brien et al. 2006; WiUingale et al. 2007; Liang et al. 2007). This shallow 
decay phase is obviously beyond understanding of the standard afterglow model (Meszaros 
& Rees 1997; Sari et al. 1998), but is generally accepted to be due to continuous energy 
injection into relativistic blast waves (e.g., Zhang et al. 2006; Nousek et al. 2006; Fan & 
Xu 2006; SoUerman et al. 2007; de Pasquale et al. 2006, 2007). In particular, Liang et al. 
(2007) recently compared the closure relations derived from a simple energy- inject ion form, 
i^inj oc t""^, with the observed temporal and spectral indices of afterglows of 53 GRBs and 
then argued that a roughly constant injection luminosity may be favored by the observations. 
However, although this simple injection form is usually assumed in model fitting, the specific 
nature of the energy injection should be concerned and studied further. 

The physical nature of an injecting energy flow discussed in the literature involves two 
possible candidates: (i) ejecta that consists of a series of shells with a distribution of Lorentz 
factors, which are dominated in energy by baryons (Rees & Meszaros 1998; Panaitescu et al. 
1998; Sari & Meszaros 2000; Granot & Kumar 2006; Liu et al. 2007), and (ii) a post-burst 
energy flow that results from a long-lasting activity of the central engine (Dai & Lu 1998a, 
1998b; Zhang & Meszaros 2001a; Wang & Dai 2001; Dai 2004; Fan & Xu 2006; Yu & Dai 
2007). In the former case, the total energy carried by the ejecta is released from the central 
engine during the prompt phase of a GRB. Subsequently, during the afterglow phase, lower- 
velocity shells catch up and collide with foregoing, higher-velocity but decelerated shells 
continuously because of an assumed power-law distribution of Lorentz factors in the ejecta. 
This persistent collision leads to a long-lasting Newtonian or trans-relativistic reverse shock 
that propagates into the ejecta. Meanwhile, the shock-heated ejecta pushes the outward- 
moving relativistic blast wave and thus reduces the deceleration of the blast wave effectively, 
which accounts for the shallow decay. In the latter case, most of the energy that produces 
an afterglow is argued to be continuously released after (not during) the burst, but the GRB 
ejecta may provide only a relatively small amount of energy for the blast wave. Initially, the 
post-burst energy flow may be dominated in energy by the Poynting flux. However, as it 
propagates outward, the flow would evolve to an ultrarelativistic kinetic-energy flow domi- 
nated by a component of electron-positron pairs through some mechanisms (e.g., magnetic 
reconnection) at larger radii (Coroniti 1990; Michel 1994; Kirk & Skjaeraasen 2003). It is 
this ultrarelativistic leptonic wind (rather than the electromagnetic flux) that feeds the blast 
wave and maintains a long-lasting relativistic reverse shock. The emission of the more and 
more energetic blast wave and the emission of the reverse shock together give rise to the shal- 
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low decay (Dai 2004). Therefore, we conclude from the above discussion that the injecting 
energy flows are likely to be matter-dominated. Furthermore, based on the component that 
is dominative in energy, the injecting flows perhaps have two types: a baryon-dominated 
outflow and a lepton-dominated outflow. The two corresponding representative models as 
described above are here called the radially structured ejecta (RSE) model (Rees & Meszaros 
1998) and the relativistic wind bubble (RWB) model (Dai 2004), respectively. 

In order to provide an effective test on the models mentioned above, a careful inves- 
tigation of high-energy emission is important and urgent, as the launch and detection of 
GLAST are upcoming (Ritz 2007). In this paper, we calculate high-energy afterglow emis- 
sion during the shallow decay phase of X-ray afterglows based on these two models and 
give corresponding observational signatures. Recently, Wei & Fan (2007), Fan et al. (2007), 
and Gou & Meszaros (2007) have made some attempts on high-energy afterglows. In their 
papers, the authors studied only the effect of the forward-shocked medium with the simple 
energy-injection form of i?inj oc t~'^, but didn't consider the effect of the injecting flow itself. 
However, as found by Dai (2004), Uhm & Beloborodov (2007), Genet et al. (2007), Yu & Dai 
(2007), and Liu et al. (2007), the reverse shock that propagates into the injecting flow could 
play an important role in the emission in X-ray and/or optical bands in some situations. 
Therefore, high-energy emission features due to the reverse shock are also expected. In the 
RSE and RWB models, the high-energy photon emission from shocked materials is mainly 
produced by inverse Compton (IC) scattering of shock-accelerated electrons (for simplicity, 
electrons and positrons are not differentiated) off synchrotron seed photons (Sari & Esin 
2001). Besides the synchrotron self-Compton (SSC) radiation from two shocked regions, 
we also consider two combined IC (GIG) processes, i.e., scattering of reverse shock photons 
by forward-shocked electrons and scattering of forward shock photons by reverse-shocked 
electrons, as pointed out by Wang, Dai & Lu (2001). 

The structure of this paper is organized as follows: in the next section we provide a 
unified description for the dynamics of fireballs in the RSE and RWB models. In §3, we 
present the energy distributions of shocked electrons that are determined by the shock- 
acceleration effect and the synchrotron & IG cooling effect, and formulate calculations of the 
synchrotron and IG radiation (including SSG and GIG). In §4, we show numerical results 
of the dynamics, spectra and light curves for some typical parameters in the two models. 
Finally, in §5, a summary is given and the observability of high-energy emission by the Large 
Area Telescope (LAT) instrument of GLAST is discussed. 
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2. Dynamics 

As illustrated in Figure 1, in both the RWB and RSE models, the system can be divided 
into four regions by the forward and reverse shock^: (1) the unshocked ambient medium, (2) 
the forward-shocked medium, (3) the reverse-shocked materials (i.e., shocked leptonic wind 
for the RWB model or shocked GRB ejecta for the RSE model), and (4) the unshocked cold 
wind or GRB ejecta, where regions 2 and 3 are separated by a contact-discontinuity surface. 



2.1. Structure of injecting flows 

Figure 1 also shows illustrative Fj-distributions in all the regions. The main differences 
between the two models are in the structure and composition of injecting flows: 

(i) The RWB model. Following the analysis of Dai (2004) for a magnetar-driven wind, 
we simply assume that the leptonic wind propagates outward with a constant luminosity L„ 
and a constant bulk Lorentz factor F^ during a period of after the burslEI. Thus, we can 
calculate the number density [n'f^^^) and the bulk Lorentz factor (F^^^^^) of the pre-shock 
materials at the reverse shock front as 

/(RWB) _t(T3\_ -^w ,^ X 

= F., (2) 

where me is the rest mass of electrons and R is the radius of the system in the thin shell 
approximation. 

(ii) The RSE model. As suggested by Rees & Meszaros (1998), the mass distribution in 



^It takes only a short time (tens to hundreds of seconds) for a reverse shock to cross the ejecta of a 
typical GRB in the RWB model. Meanwhile, two forward shocks forming initially during the interaction of 
the ejecta both with the medium and with the leptonic wind are assumed to merge to one forward shock. In 
addition, the contribution of the GRB ejecta to the afterglow emission should be negligible as compared with 
the shocked medium and the shocked wind during the emission period of our interest. Thus, for simplicity, 
the structure of an RWB can be described approximately by Figure 1(a). However, for the RSE model, the 
situation is different. The reverse shock propagates into the GRB ejecta for a long time and thus influences 
the afterglow emission during the first hours significantly (see numerical calculations in §4). 

quantity Q' with a prime is measured in the comoving frame of a certain region, while the other 
quantities Qs are measured in the local medium's rest frame except for specific declarations. In addition, a 
subscript "i — 1,2,3,4" represents a quantity describing region "z". 
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the wide GRB ejecta associated with a distribution of bulk Lorentz factors reads 

Mej(>rej)«r-. (3) 

Moreover, we assume that the Fej-distribution in the GRB ejecta satisfies 

rej(x)0CX-^^ (4) 

where x is the displacement of the reverse shock propagating into the ejecta. If the Lorentz 
factor of the end of the ejecta (i.e., the minimum bulk Lorentz factor) is denoted by Fej^min, 
then the total displacement of the reverse shock when it crosses the ejecta could be estimated 
as /?cross/2rgj jj^ij^ ~ 10^'^ — 10^^ cm. Under these assumptions, the associated injecting energy 
distributes with respect to x as E^^ oc x*/'', and n'^^^^^-* and r4^^^^ are given by 



./(RSE) _ ^/ _ dM^j/dx 

47ri?2rej(x)mp' 



n-, ' = n',^{x) = , (5) 



rr^) = rej(x), (6) 

where nip is the proton rest mass. 



2.2. Dynamic equations 

Now, we describe the dynamic evolution of the systems under the effect of the two types 
of injecting flow. A resultant long-lasting reverse shock transforms the injecting energy into 
the internal energy of the reverse-shocked materials continuously. Meanwhile, a part of fresh 
energy is further transformed to the kinetic energy (£'k,2) of region 2 through the work done 
by region 3: 5W — AixP^R^dR, where is the pressure of region 3. Then, we have 

8W - dEy,^2 = d [{Tl - l)M3wc2] , (7) 

where M^^ is the rest mass of the swept-up medium and T2 is the average bulk Lorentz 
factor of the shocked medium. Thus, the dynamic evolution of region 2 is described by (Dai 
2004) 

^ ^ 47ri?^ [PUc^ - {Tl - l)nimp] 

dR IT^M^^ ' ^ 

where dMg^/ dR — A^nR^nirrip is applied and rii is the proton number density of the ambient 
medium. To integrate the above equation, the pressure of region 3 should be calculated by 

^3 = ^(r34 - l)(4r^4 + 3)r^XeC^ (9) 
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where = r3r4(l — PsPa) is the Lorentz factor of region 3 measured in region 4 and (3i 
is the velocity in units of c, and rrire represents the electron and proton rest masses for the 
RWB and RSE models, respectively. The evolution of can be obtained from the equality 
of Lorentz factors of the two sides of the contact discontinuity surface, assuming that the 
shocked medium in region 2 satisfies the adiabatic self-similar solution of Blandford & McKee 
(1976) 

Ts = T,x~'^'. (10) 
To fix the self-similar variable we use the following relationship, 

= ^Tlmm,c\-''/'\ (11) 

Combining Eqs. [HI [101 and[TT]to eliminate x? we can solve and as functions of r2 and 
thus integrate Eq. [HI Furthermore, the rest mass of region 3 is obtained from 

^ = A7TR\f3, - /3Rs)r4n;mre, (12) 

where P^s = (Xs^sPs ~ ^in'^^Pi) / (Xs^'s ~ ^4^'^) is the velocity of the reverse shock. 

When the reverse shock crosses the wind or GRB ejecta, the dynamic evolution of 
region 3 is described by Kobayashi & Sari (2000) and Kobayashi (2000). Then, the pressure 
of region 3 in Eq. [H] decreases significantly and thus is negligible. As a result, the dynamic 
equation of region 2 returns to the form of the standard afterglow model in Huang et al. 
(1999). 



3. Electron distributions and emission mechanisms 



After the dynamic evolution equations are given, the internal physics of the shocked 
regions in the RWB and RSE models can be considered as follows. 



3.1. Electron energy distributions 

As the forward and reverse shocks propagate, the bulk kinetic energies of the shells 
are gradually transformed into the internal energy of the shocked materials, the density of 
which is denoted by e^. The internal energy will be partially carried by the accelerated 
electrons and magnetic fields, the energy densities of which are factions ee,i and ee.i of 
e-, respectively. Through the shock acceleration, the electrons behind the forward/reverse 
shock will obtain an initial energy distribution A^/niilT't) with the minimum Lorentz 
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max, I ^ 



factor = Xi[{p — 2)/{p — l)]ee,j(rj — 1) and the maximum Lorentz factor 7^ 

10®[i?-(l + Fj)]"^/^, where A2 = m^/m^, A3 = m^c/^e, "^2 = ^2, fs = T'^^, and B'- is the 
magnetic field strength. The occurrence of the Compton parameter Yi is induced by the IC 
coohng effect. By considering the coohng effect of the synchrotron and IC radiation, we can 
obtain the actual electron energy distributions (Huang et al. 2000): 
for 7'^_i < 7'^_i, 



at' f^J \ ) ^ i if T c,i — T i — T m,i; /-i q\ 

y-p-1 if Y .<y.<y .- ^^3) 

II '■'■I m,i ^11 — / max, 11 



and for 7'^ • > 7' 



m,i' 



max,« • 

The cooling Lorentz factor 7'^, ^ is defined by equaling the cooling timescale and dynamic 
expansion timescale of the system, and it reads 

, _ 67rmcC 

with (Tt being the Thomson scattering cross section. Here, the dynamic timescale measured 
in the observer's frame in calculating 7'^, is assumed to be approximately equal to the ob- 
server's time t since the trigger. In addition, the self-absorption of the synchrotron photons 
is ignored in Eqs. [13] and [T3] for the high-energy emission that we are interested in. Follow- 
ing Sari & Esin (2001), the Compton parameter Yi in Eq. [151 defined as the ratio of the IC 
(including SSC and CIC) luminosity to the synchrotron luminosity of electrons, is estimated 
by 



Tl 

_ ^IC,2 


TJi _L Ujl 
'^syn,2 ^ 2 syn,3 


~ Tl 

syn,2 


K,2 


Tl 


Tf _L iTJi 
•^syn.S ^ 2 syn,2 


-^syn.S 





V — i<^,i sjn,A ■ 2 syn,z „n 
^3 = 77 - 77, , \^') 

where Ul.^^ ^ = ?7jee,ie^/(l + Yi) and U^ ^ = eB.iC- are the energy densities of synchrotron seed 
photons and magnetic fields, respectively. The radiation efficiency rji reads 



1> for I'ci < 7 'm,i, 

(7 c,i/ 7 m,j) ^5 for 7 c,j ^ 7 m,r 



A factor of | in Eqs. [16] and [T7| occurs, because only about one-half of seed photons from one 
shocked region will diffuse into the other one for the CIC process. In two extreme situations. 
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equations [16] and [T7| can be simplified as 



i) for f/;,_3 > f/,;^^2, 



Y, = rjsp^{l + Y,)-\ Y, = r^s^{l + Ysr'; (19) 
^^6,262 eB,3 



ii) for [/^'y^ g < f/^'y^ a, 



Y, = r^,^il + Y,)-\ Y, = V2f^il + Y,)-\ (20) 



3.2. SSC and CIC emission 

Once tlie electron distribution is known, tlie synclirotron emissivity of electrons in region 
i at frequency z/' is calculated directly by (Rybicki & Liglitman 1979) 



= I ''^-'''VaV.)^ ( - I . (21) 



where qe is the electron charge, u'^ = dt'^'lq^B'J {ATrm^c) , ^(u) = u K^jj,{k)dk, and K^j^^ik) 
is the Bessel function. 

Accompanying with the synchrotron radiation, the electrons also lose their energy 
through upscattering the synchrotron seed photons. As usual, the first-order IC scatter- 
ing is considered and the higher-order processes are neglected. According to Blumenthal 
& Gould (1970), when the Klein-Nishina suppression is considered, the IC emissivity (at a 
frequency v') of electrons in region i upscattering seed photons from region j is calculated 
by {i = j for SSC and i j for CIC) 

ef^^\u') = 3aT / diMl'^) [ du]^^g{x,y), (22) 

where g{x,y) = 2y\ny + {I + 2y){l - y) + |f^(l - y) with x = Ay^hu] / rrieC^ , y = 
hv' /[x{^[mf,c^ — hv')]. z/j and /| are the frequency and the corresponding flux density of 
the incident photons from region j, respectively, which are measured in the comoving frame 
of region i. 

The observed synchrotron and IC flux densities at a frequency v (measured in the 
observer's frame) from region i are given respectively by (Huang et al. 2000) 

F7m= Fdevf^'-T^^ (23) 
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sine sT'^'^Viu) 
2 A-nDlV^ 



(24) 



where is the luminosity distance and Vi = Ti{l — f3iCos6) is the Doppler factor. Here, 
the viewing angle is assumed to be zero and 6 is the angle between the moving direction of 
the emitting material and the line of sight. The volumes of the shocked regions in Eqs. 
[23] and [2l] are functions of 6 because the integrations are performed over the equal-arrival 
surface. In addition, the jet correction should be considered if the GRB ejecta is coUimated. 



For the RWB and RSE models, we first calculate the dynamic evolution with radius. 
The initial value of R is taken to be ~ 10^^ cm as the deceleration radius. For the RWB 
model, as in Yu & Dai (2007), the isotropic wind luminosity is = 4 x 10'^'^ Bf^ erg s~\ 
the wind duration = 5 x 10'^ B^^ s, and the bulk Lorentz factor of the wind = 10^, 
where B^ is the magnetic field strengthen of the central magnetar in unit of 10^^ G. Here, 
we take B14 = 4 for a typical magnetar. Thus the total energy of the wind is = 2.0 x 10^^ 
erg. However, the isotropic kinetic energy of the previous GRB ejecta is relatively small, 
^RWB _ ]^g5i gj,g_ Pqj, consistency, we assume that, in the RSE model, the isotropic kinetic 
energy carried by the wide GRB ejecta is E^^^ = 2.1 x 10^^ erg, and the distribution-related 
parameters are Fej,min = 30, Fej,max = 500, and s = b = 1.5. In addition, the number density 
of the ambient interstellar medium is ni = 1 cm~^ in both cases. 

The dotted lines in Figure 2 represent the dynamic evolution of blast waves with a 
certain energy (i.e., without energy injection), and the solid lines show the energy injection 
case. We can see clearly from this figure that the energy (oc F2) of the shocked medium 
increases gradually until the energy injection is over. As analyzed in Yu & Dai (2007), the 
final energy carried by the shocked medium is a fraction ~ 67% of the total injecting energy 
for the RWB model versus ~ 90% for the RSE model. Meanwhile, the other fraction of the 
injecting energy (~ 33% and ~ 10% for the RWB and RSE models, respectively) should be 
shared by the reverse-shocked material. In addition, we can know that the reverse shock 
is relativistic in the RWB model, but Newtonian or trans-relativistic in the RSE model, 
according to an estimation of the Lorentz factor of region 3 relative to region 4 from 



with the values of F3 as shown in Figure 2. This difference between the reverse shocks in 



4. 



Numerical Results 



4.1. Dynamic evolution 




(25) 
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the two models is just the reason why the reverse shock in the RWB model can share more 
injecting energy than the one in the RSE model. 

4.2. Spectra and light curves 

To calculate the emission of shocked materials, we take the microphysical parameters 
p — 2.3 and es = 0.01 for all shocked regions. Then, is calculated by ^/e^ for the baryon- 
dominated regions (i.e., regions 2 & 3 in the RSE model and region 2 in the RWB model) as 
argued by Medvedev (2006) and (1 — ee) for the lepton-dominated region (i.e., region 3 in the 
RWB model). The luminosity distance of a GRB is taken to be D^, — 1 Gpc, corresponding 
to a redshift of 0.2. 

Figure 3 shows the hght curves of an RWB. From this figure, we obtain the following 
results. First, the shallow decay phase of an X-ray (kcV) afterglow is produced during the 
early hours (left upper panel), which is mildly dominated by the reverse shock emission 
(for a detailed discussion see Yu & Dai 2007). The synchrotron radiation is the dominative 
mechanism in this band. Second, an obvious hump accompanying with the X-ray plateau 
occurs in a high-energy (MeV & GeV) afterglow, which is dominated by the IC (especially 
SSC) emission from the reverse-shocked wind. Finally, a comparison of light curves in 
different bands (right lower panel) shows that the high-energy (especially GeV) emission 
flux is mildly or even significantly higher than the one in X-ray band. This feature is also 
implied obviously by the spectra as shown in Figure 4. The peak flux of the IC emission at 
~ 1 GeV is about one order of magnitude higher than the synchrotron peak at ~ 10 eV. 
Moreover, the IC component is dominative above sub-McV and, in all bands, the emission 
from the reverse-shocked wind is more important than the one from the forward-shocked 
medium at 1000 s. 

The results for the RSE model are shown in Figures 5 and 6. From these two figures, we 
find that, first, the X-ray shallow decay phase during the early hours is also produced in this 
model as shown in left upper panel of Figure 5, which is dominated by the forward shock. 
Second, the high-energy afterglow light curves also have a flattening segment. However, 
differing from the RWB model, the contribution to the high-energy plateau from the IC 
emission plays a role only in GeV band, whereas the MeV emission is totally contributed 
by the synchrotron mechanism during the total afterglow phase. This result arises from the 
relative weakness of the IC component in the RSE model as shown in Figure 6. Therefore, 
the flux in MeV or GeV band is lower than or approximately equal to the X-ray flux as 
shown in right lower panel of Figure 5. In addition. Figure 6 shows that the contribution 
of the emission from the reverse-shocked ejecta is negligible in all bands for the parameters 
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that we adopt. 

In conclusion, we obtain very different high-energy afterglows in the RWB and RSE 
models with a same amount of injecting energy that gives rise to similar X-ray afterglows. 
The RWB model predicts more significant (about one order of magnitude stronger) high- 
energy emission than the RSE model. The reasons for this difference are that the reverse 
shock in the RWB model is relativistic and that the energy of the reverse shock is almost 
totally carried by electrons. The former reason enables the reverse shock to share more (i.e., 
33% vs 10%) injecting energy as discussed in Sect. 4.1, and the latter reason dramatically 
increases the emission efficiency of the reverse-shocked material through the enhanced IC 
emission. Therefore, a higher fraction of injecting energy can be radiated from the shocked 
regions in the RWB model. In the baryon-dominated injection model, however, this sig- 
nificant enhanced IC component doesn't happen even if the emission fiux from the reverse 
shock exceeds the one from the forward shock under some extreme conditions, because most 
of the injecting energy in both shocked regions is locked in the baryons whose emission is 
weak. Therefore, we argue that the difference in high-energy emission between the RWB and 
RSE models is essentially due to the physical distinction between the two types of energy 
injection. Thus, reasonable variations of the parameters for the models should not change 
our results significantly. 



5. Summary and discussion 

The discovered shallow decay phase of GRB X-ray afterglows suggests long-lasting en- 
ergy injection into relativistic blast waves. The injecting flow is likely to be dominated in 
energy by the component of either baryons or leptons. In this paper, we first provide a uni- 
fied description for dynamics and radiation in two representative models, i.e., the RSE and 
RWB models. Through maintaining a long-lasting reverse shock and doing work to forward- 
shocked medium, both types of energy injection can produce the flattening segment in X-ray 
afterglow light curves easily. Second, we pay attention to calculations of the simultaneous 
high-energy emission that is due to synchrotron and, especially, IC (including SSC and CIC 
processes) radiation from the shocked materials. Our results show that, during the shallow 
decay phase of X-ray afterglows, there is a plateau (even a hump) in high-energy light curves 
in both the RSE and RWB models. As argued by Wei & Fan (2007), we suggest that the 
plateau/hump might account for the delayed high-energy emission of some bursts such as 
GRB 940217. 

We also find that the high-energy emission derived from the two models has different 
observable features, e.g., morphologies of the light curves and spectra. In particular, more 
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significant high-energy emission is predicted by the RWB model, because more injecting 
energy in this model is carried by electrons and thus the IC scattering is enhanced signifi- 
cantly. This difference in high-energy emission between the two models could be tested in 
the upcoming GLAST era. As proposed by Zhang & Mcszaros (2001b) and Gou & Meszaros 
(2007), the fluence threshold of the LAT instrument aboard GLAST for an effective observa- 
tion time tes in an observation energy band centered around a photon energy E is estimated 
by -Fthr ~ ['^'o(^/^efr)^''^]-E'teff, where the flux sensitivity, [$o(2"/^eff)^'^^], scales as due 
to a limitation by the background for a long-time observation and $o is the sensitivity over 
an exposure time T. Specifically, over a one- year exposure period, the sensitivity above 100 
MeV is ~ 1.33 x 10~^ ph s~^ cm~^ for point-source observations of GRB afterglows (Gou & 
Meszaros 2007). On the other hand, for a short-time observation, the fluence threshold is in- 
stead calculated by Fthr = ^E/Acs with the assumption that at least 5 photons are collected 
on an effective area A^s of the instrument. Taking E = 400 MeV and A^q = 6000 cm~^, we 
obtain the fluence threshold of GLAST LAT, 

_ f 5.3x10"^ ergcm-^ f < 2.43 x lO^s, 

~ \ 3.4 X 10-9^1/2 erg cm-^, t > 2 AS x lOS, ^ ' 

where an observation efficiency r) = tes/t — 0.5 due to the occultation by the earth is 
considered. With this fluence threshold, the observabihty of high-energy emission in the 
RWB and RSE models by GLAST LAT is shown in Figure 7. It can be seen that both 
the high-energy emissions predicted by the two models can be detected by GLAST LAT for 

the given parameters. The partial fluence of the high-energy emission in the RWB model 
is about an order of magnitude higher than that in the RSE model for a same amount of 
injecting energy. Here, the partial fluence is defined as an integration of the fiux density 
(F^) over the GLAST LAT energy band [20 MeV, 300 GeV] and the time intervals [0.5t, t] 
as done in Gou & Mcszaros (2007). 

Although the detection of the RWB model by GLAST LAT is obviously easier, joint 
observations of GLAST and Swift are necessary to distinguish between the RWB and RSE 
models. Simultaneous observations of a high-energy afterglow by GLAST LAT and an X- 
ray afterglow by Swift XRT for a same GRB could provide a ratio of fluxes (z^iv) in the 
high-energy band to the X-ray. As imphed by the right lower panels of Figures 3 and 
5, for example, the injecting flow could be dominated in energy by baryons if the ratio 
{i^F,^)QeV / {^F^)keV is closc to 1. On the other hand, if the ratio is signiflcantly higher than 
1, there could be lepton-dominated energy injection. 
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Reverse-Shocked Forward-shocked Reverse-Shocked Forward-shocked 

leptonicwind medium GRB ejecta medium 




r r 
(a) RWB model (b) RSE model 



Fig. 1. — Illustrations of four regions divided by the forward and reverse shocks in the RWB 
(left) and RSE (right) models, and the corresponding distributions of Lorentz factors (not 
scaled) . 
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Fig. 2. — Evolution of the Lorentz factors of regions 2 and 3 as the shells in the RWB (left) 
and RSE (right) models expand. A dotted line represents the dynamic evolution of the shell 
with a fixed amount of energy as labeled. 
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Fig. 3. — Light curves in X-ray (keV) and high energy gamma-ray (MeV, GeV) bands for 
the RWB model The emission from the forward-shocked medium is represented by thick 
hues, while the reverse shock corresponds to thin hues. The sohd, dash-dotted, and dashed 
hues are due to SSC, CIC, and synchrotron processes, respectively. 
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Fig. 4. — Synchrotron and IC spectra at 1000 s for the RWB model. The sohd hnes represent 
the emission contributed by the forward-shocked medium, while the dotted lines correspond 
to the reverse-shocked wind. 
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Fig. 5.~- 



Same as in Figure 3 but for the RSE model. 
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Fig. 6. — Same as in Figure 4 but for the RSE model. 
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Fig. 7. — The partial fluence (defined as J^'^ F^dvdt, where hui — 20 MeV and hi'2 — 300 
GeV) curves for the RWB and RSE models. The dotted line represents the fluence threshold 
of the GLAST LAT at 400 MeV. 



